Rigidity dependent knee and cosmic ray induced 



Julian Candia a and Esteban Roulet & 



high energy neutrino fluxes 

o ' 
a, 

^ '. a Departamento de Fisica, Universidad Nacional de La Plata, CC67, 

La Plata 1900, Argentina 
b CONICET, Centro Atomico Banloche, Av. Bustillo 9500, 
Bariloche 8^00, Argentina 



February 2, 2008 



I 

i-^ ■ Scenarios in which the knee of the cosmic ray spectrum depends on the par- 

^p" 1 . ticle rigidities usually predict that the cosmic ray composition becomes heav- 

ier above the knee and have associated a change in the spectral slope of each 
individual nuclear component which is steeper than the change (Aa ~ 0.3) 
observed in the total spectrum. We show that this implies that the very high 
energy (E v > 10 eV) diffuse neutrino fluxes produced by cosmic rays hitting 
the atmosphere or colliding with the interstellar medium in the Galaxy will be 
significantly suppressed, making their detection harder but also reducing the 



background for the search of other (more challenging) astrophysical neutrino 
sources. 



1 Introduction 

The study of atmospheric neutrinos of sub-GeV and multi-GeV energies has been of 
paramount importance in the recent past, providing in particular the first clear evi- 
dence in favor of neutrino oscillations, and hence of non- vanishing neutrino masses. 
These neutrinos are produced mainly in the decay of pions and kaons (e.g. 7r~ — > 
fi~V^, with ji" — > e~V e Vy) produced by the cosmic rays (CRs) hitting the upper 
atmosphere and generating air showers with a significant hadronic component. Due 
to the interplay between the mesons' decay and their interactions in the air, as 
the energy increases (with the associated relativistic dilation of their decay times) 
the mesons lose an increasing amount of energy before they decay, and also the 
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muons produced can reach the ground before decaying. As a consequence of this, 
the predicted atmospheric neutrino spectrum results steeper than the CR spec- 
trum typically by an extra power of the energy, with dN u /dE u oc E~ 3 - 7 , while 
dN CR /dE CR oc Ec 2 R 7 [1, 2, 3]. 

Comparing the charged mesons' decay length, L = 7cr, where in particular 
L n ~ 5.6 km(E n /100 GeV) and — 7.5 km(i?fc;/TeV), it is clear that above 
100 GeV the pions traverse several attenuation lengths (corresponding to approxi- 
mately 120 g/cm 2 ) in the atmosphere before decaying, and due to this the neutrino 
fluxes above a few hundred GeVs actually arise mainly from K decays. Kaons in 
turn are also significantly attenuated before decaying for energies above the TeV, 
with the effect that the prompt neutrinos from charmed particle decays become 
increasingly important. 

Several groups have studied the charm particle production by CRs [4, 5, 6], which 
actually requires to take into account next to leading order (NLO) processes, which 
increase the charm production by more than a factor of two, and the cross sections 
are sensitive to the values of the partonic distribution functions at very small values 
of the scaling variable x, what introduces further uncertainties in the results, due to 
the need to extrapolate beyond the measured range. Even if the charmed mesons 
(D ± , D°, D s , A c , ...) are produced at a rate ~ 10~ 4 with respect to the non-charmed 
ones, their very short decay time, with Lp ~ 2 km(.Ex>/10 PeV), implies that the 
prompt neutrino spectral slope just follows the CR slope, and hence the prompt 
neutrinos eventually dominate the atmospheric neutrino fluxes above a PeV in the 
vertical direction, and somewhat above that energy at large zenith angles 1 . 

It is important at this point to take into account the fact that the primary 
CR spectrum shows a steepening at the so-called knee, corresponding to Ek nee ~ 
3 x 10 15 eV, with dNcn/dEcR oc Eq R for Eqr > Eknee- The associated steepening 
in the induced neutrino fluxes has been obtained in the literature [4, 7], but these 
works usually assume for simplicity that the dominant CR component consists of 
protons, and only some crude attempts have been done trying to generalise these 
predictions to heavier compositions [6]. 

Another contribution to the diffuse neutrino fluxes produced by CRs are those 
resulting from the interactions of the CRs present all over the Galaxy with the 
ambient gas in the inter stellar medium (ISM) [8, 9, 10, 7]. The very low densities 
in the ISM (n IS M — 1/cm 3 ) imply that essentially all mesons produced decay in this 
case without suffering any attenuation, and hence the neutrino fluxes are just the 
conventional ones from ir, K and [i decays, and moreover their spectral slopes follow 
the behaviour of the primary CR spectrum. These neutrino fluxes are proportional 
to the column density of the ISM along the particular direction considered, and 
hence the flux is maximal in the direction to the Galactic Center, it is generally 

1 At large zenith angles the mesons produced high in the atmosphere are less affected by the 
attenuation before decay due to the fact that the atmosphere they traverse is more tenuous. 
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enhanced near the galactic plane while it is minimum in the direction orthogonal to 
it. 

Along directions near the plane of the Milky Way, these diffuse fluxes of galac- 
tic neutrinos may overcome the atmospheric neutrino background at energies larger 
than ~ 10 14 eV (the actual energy depends on the model adopted for charm pro- 
duction and on the assumed ISM column density, see below), but in the direction 
orthogonal to the plane they remain below the expected flux from prompt atmo- 
spheric neutrinos from charmed particle decays up to the highest energies [7]. Here 
again the usual simplifying working hypothesis is that the CR flux consists mainly of 
protons, and only crude attempts have been done trying to generalise the predictions 
to heavier compositions [9]. Although some works have found a strong sensitivity to 
the assumed composition both for atmospheric and galactic CR induced neutrinos, 
this issue is generally not treated appropriately and hence a more thorough anal- 
ysis is required in order to estimate the neutrino fluxes with some confidence for 
E v > 10 14 eV. 

The very high energy neutrino fluxes produced by CRs both in the atmosphere 
and in the Galaxy are one of the important targets of present and future neutrino 
telescopes, such as AMANDA, BAIKAL, ANTARES, ICECUBE, etc.. They are 
interesting per se, but they also represent the main background in the search of 
the fluxes arising from other potential neutrino sources, such as Active Galactic 
Nuclei (AGNs) [11] and Gamma Ray Bursts (GRBs) [12], which might give rise to 
detectable fluxes just in the energy range 10 14 -10 16 eV [13]. It is hence of primary 
importance to establish with some confidence the actual value of the CR induced 
neutrino fluxes, and the aim of the present work is to point out that the predictions 
are quite sensitive to the detailed composition and to the behaviour of the spec- 
trum of the individual nuclear components of the CRs with energies above those 
corresponding to the knee in the spectrum. 

A major problem that one has to face here is that although the existence of the 
knee has been known for more than forty years, there is still no consensus on the 
underlying physics responsible for this feature, and different proposed explanations 
for it can lead to quite different predictions for the behaviour of the CR composition 
and the individual spectra. On the observational side, due to the indirect nature 
of the measurements, which are based on the analysis of extensive air showers, also 
the situation is far from settled. 

A large class of possible scenarios to explain the knee are based on a rigidity 
dependent effect, such as a change on the CR acceleration efficiency at the sources 
[14, 15, 16] or on a change in the escape mechanism of CRs from the Galaxy [17, 18, 
19, 20], both of which, being magnetic effects, depend on the ratio Ecr/Z, where 
Z is the CR charge. Hence, in these scenarios the different nuclear components 
essentially steepen their spectra at an energy E z — ZE knee , and in this way the light 
components become more suppressed at smaller energies than the heavy ones. This 
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leads to the prediction that the CR composition should become increasingly heavier 
above Eknee, in agreement with the latest observations of KASCADE [21] and the 
EAS-TOP/MACRO experiments [22]. In addition, the change in the spectral slope 
of each individual component usually turns out to be steeper than the change in 
slope of the total observed spectrum, which is Aa ~ 3 — 2.7 = 0.3. In particular, in 
the so-called diffusion/drift model [19, 20, 23, 24] in which the knee is due to a change 
in the escape mechanism of CRs from the Galaxy from one dominated by normal 
diffusion (with diffusion coefficient D oc E 1 ^ 3 ) to one dominated by drift effects 
(with D oc E), the spectral slope of each individual nuclear component changes by 
an amount Aa z — 2/3. 

As we will show below, these rigidity dependent scenarios have a profound im- 
pact on the predictions of high energy diffuse neutrino fluxes, and hence a clear 
understanding of the physics responsible for the knee and also better measurements 
of air showers at these energies are crucial to predict the expected neutrino fluxes 
and to extract conclusions from them once they will be measured. 



2 The cosmic ray spectrum 

Figure 1 shows the observed cosmic ray differential spectrum from several experi- 
ments [25]. As mentioned above, the total CR spectrum is well described by power 
laws dNcn/dEcR oc E^, with the spectral index changing from a ~ 2.7 to a ~ 3 
at the so-called knee, occurring at E knee ~ 3 x 10 15 eV. The CR spectrum turns 
harder again at the ankle, corresponding to E ank ie ~ 5 x 10 18 eV, a feature usually 
explained as the crossover between the dominance of the galactic CR component 
below the ankle and of the extragalactic one above it. 

A simple way of describing the galactic CR flux within a rigidity dependent 
scenario is to assume that each galactic component of charge Z is given by 

= <?z - 



dE ™ 0< + 0>' 
where (f> z (0^) is the CR flux below (above) the knee, and is given by 



!>< I „ , fE\- az j i 



In this expression, O = 3.5 x 10 _13 m _2 s _1 sr _1 eV _1 is the total CR flux at the 
reference energy E , hereafter adopted as E = 1 TeV, f z is the fractional CR 
abundance at the same energy, a z is the low energy spectral index, Aa is the spectral 
index change across the knee (which is the same for all CR components) and E k is 
a parameter that fixes the position of the knee, which is conveniently chosen to fit 
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Table 1: Cosmic ray fractional abundances (for E = 1 TeV) and low energy spectral 
indices from hydrogen to nickel. 

the observations. The expression in Eq. (1) provides then a smooth interpolation 2 
between the flux 0^ at low energies (E < ZE k ) and 0^ at high energies (E > ZE k ). 
The CR abundances at E — 1 TeV and the low-energy spectral indices are shown 
in Table 1, and were taken from the data compiled in [26, 27]. 

For the extragalactic component, we will assume for definiteness that it consists 
of protons (this would anyhow only affect the predictions for E u > 10 17 eV) and 
that it is given by 

0,r 6 ' 8xir34 (l^v)' 2 ' 4m " V ' sr " leV " 1 ' (3) 

i.e. a flux similar to that considered in [24]. 

For the purpose of illustration, we display in Figure 1 the cases of both Aa = 2/3, 
corresponding to the diffusion/drift model (when a Kolmogorov spectrum of fluctu- 
ations is assumed for the random magnetic field component in which CR particles 
propagate [20, 23, 24]), and Aa = 2, corresponding to an extreme model advocated 
in [27]. Figure 1 shows the fits that result from taking Aa = 2/3 with the parameter 
value E k = 2.2 x 10 15 eV, and Aa = 2 with E k = 3.1 x 10 15 eV. While the former 
exhibits an excellent agreement with the experimental observations, it can be seen 
that the latter shows an abrupt suppression above ~ 10 17 eV which results from the 
steep suppression of the different CR components above their knees [27], and due to 
this we will hereafter just focus on the Aa = 2/3 case. 

2 One may eventually generalise the expression in Eq. (1) as <j> = (f> K • 4> > /{4> <n + <j) >n ) 1 ^ n , where 
the parameter n will control the 'width' of the knee, with n»l leading to an abrupt change in 
slope. We find however that n = 1 provides a very satisfactory fit to the CR data for the Aa = 2/3 
case considered. 
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Figure 1: CR spectra obtained from considering rigidity-dependent scenarios 
parametrized according to equations (1) and (2). For the spectral index change 
Aa = 2/3, the main contributions to the total flux, which correspond to nuclei of 
H, He, and Fe, are also indicated. The dotted straight line corresponds to the 
extragalactic flux given by equation (3). Also shown are the relevant experimental 
observations. 

3 The flux of atmospheric neutrinos 

The CR particles reaching the top of the atmosphere produce secondary fluxes of 
hadrons and leptons that are usually described by means of coupled cascade equa- 
tions [2, 3, 4, 5], which can be solved analytically under appropriate symplifying 
assumptions. In this work we will compute the fluxes of muon neutrinos and an- 
tineutrinos, which are the ones more readily detectable at underground detectors. 
At the energies considered the neutrino flavor oscillations will have no effect on the 
atmospheric neutrinos, although they may be relevant for the neutrinos produced by 
CR interactions with the ISM (eventually redistributing in that case the produced 
u e and fluxes evenly among the three flavors, something that can be incorporated 
into the final results straightforwardly). 

The main contributions to the neutrino flux in the range of interest of this work 
(namely, for E = 10 3 — 10 8 GeV) arise, on the one hand, from the decay of the 
mesons TT ± ,K ± ,Kg and K®, which produce the so called conventional atmospheric 
neutrinos that dominate at low energies, and on the other hand, from the decay 
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of the charmed mesons D°,D ,A C and D s , which give rise to the so called prompt 
charm neutrino flux that dominates at very high energies. The muon contribution 
to the atmospheric neutrino flux is completely negligible in this context, since above 
E ~ 10 3 GeV the muons reach the ground and are stopped before decaying. 



3.1 Nucleon and meson fluxes 



jg&X) = -£(E,X) - f d (E,X) +^S kj (E,X) , (4) 



The transport equation associated to a given cascade component j can be written 

as 

\? 

where as usual X is the slant depth, Xj is the interaction length in air and \j is the 
decay length. Both characteristic lengths are measured in g/cm 2 units (i.e. they in- 
clude the factor p(X) that corresponds to the local density of the atmosphere). Also 
notice that the decay length is linear in the energy due to the Lorentz time dilation 
factor, while the interaction length exhibits instead a much milder energy depen- 
dence. The last term in equation (4) corresponds to the production/regeneration 
term given by 

S kj (E,X) = J e dE Xk(E/) — W) , (5) 

which actually couples the transport equations of different cascade components. 
Notice that a k A here refers to the total fc+air cross section, while a k A->j corresponds 
to the process /c+air — > j+ anything. Since the energies involved in the nuclear 
collisions between CRs and atmospheric nuclei are much higher than nuclear binding 
energies, only nucleon-nucleon interactions are relevant in this context. The different 
cascade components that we will take into account are hence N (that corresponds 
to the nucleon component in which protons and neutrons are grouped together) 
and M = vr ± , K ± , K° L , D°, D ± , A c and D s . As in [3], the K° s can be considered 
as contributing to the pion flux. Moreover, we will assume for simplicity that the 
only non- negligible terms will be those corresponding to regeneration (i.e. Snn and 
Smm) and to meson production by nucleons (i.e. SVm) [3, 4]. 

Let us consider in particular fluxes with the form (f) k (E,X) = E~@ k g(X) (i.e. 
we consider the separation of variables in the fluxes, in which the energy dependent 
factor has a power-law behavior). Then, it turns out that 

S k AE,X) = ^^Z kj (E), (6) 

where Z k j is the production/regeneration moment defined by 

Z k ,(E) = /"d, /.-.'^ , (7) 

Jo CT k A{E) dx 
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with x = E/E'. 

Assuming that the interaction lengths are energy independent and that the differ- 
ential production distribution is Feynman scaling, the production/regeneration mo- 
ments themselves become energy independent. Indeed, for the relevant interaction 
lengths and production/regeneration moments involving nucleons and non-charmed 
mesons, we consider the energy independent values given in [2, 3]. When necessary, 
we also take the (3^— dependence from fits to the data compiled and plotted in [2]. 
Moreover, following [5], the charm interaction lengths and regeneration moments 
are all set equal to the kaon ones. 

For the production of charmed particles, however, one has to take into account 
the energy dependent moment given by equation (7), since the results are very 
sensitive to the particular cross section adopted for the charm production process. 
Indeed, the various estimates of the prompt atmospheric fluxes calculated so far are 
found to differ by almost 2 orders of magnitude [28]. 

Due to their relatively large mass, charm quarks are usually considered to be pro- 
duced in hard processes which can be well described by perturbative QCD (pQCD). 
To leading order (LO) in the coupling constant, the processes that contribute to the 
charm production cross section are the gluon-gluon fusion process gg — > cc and the 
quark- ant iquark annihilation process qq — > cc. However, at the next to leading or- 
der (NLO) the gluon scattering process gg — >■ gg shows up, giving a very significant 
contribution to the total cross section which increases it by a factor of ~ 2 — 2.5 
[29, 30, 31]. Concerning higher order corrections, their contribution is certainly small 
since there are no qualitatively new channels opening up. Several parton distribu- 
tion functions have been proposed and provide theoretical pQCD predictions that 
fit the available accelerator data reasonably well. However, in order to calculate the 
atmospheric charm flux the parton distribution functions need to be extrapolated to 
very small parton momentum fractions, typically corresponding to x ~ 4 GeV/E [6], 
which at the highest energies is well outside of the measured region (x > 10~ 5 ), and 
the uncertainty in the extrapolation affects significantly the final results [5, 28]. In 
order to illustrate this uncertainty range, we display the results obtained with two 
different structure distribution functions, namely the CTEQ3 parton distribution 
function [32] and the Golec-Biernat, Wusthoff (GBW) model [33], which includes 
gluon saturation effects. The corresponding charm production cross sections were 
obtained either directly from the fit given in [6] (for the results based on the GBW 
saturation model), or by fitting the results given in figure 4 of [5] for the CTEQ3 
structure functions set with M = 2m c = 2/j,, with the charm mass m c = 1.3 GeV, 
and interpolating the results for different energies. For the nucleon-air cross section, 
we used the parametrization given by [34] 



a NA (E) = 280 - 8.7 In 
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Once all relevant interaction lengths and production/regeneration moments are 
determined, the nucleon and meson atmospheric fluxes can be calculated from the 
coupled cascade equations given by (4) and (6). Recalling that these expressions 
were obtained under the assumption of power-law nucleon and meson spectra, we 
will first discuss the results for an initial nucleon flux given according to 4>n(E, X = 
0) = (PonE' 1 , discussing the general case further below. 

Under the assumptions mentioned above, the nucleon flux develops indepen- 
dently from the meson fluxes and is given by 

ME,X) = e- x / AN <f )0N E-\ (9) 

where the nucleon attenuation length is defined as 

A* = j^- • (10) 

Concerning the meson cascade equations, they are usually solved by considering 
separately the low energy solution (f> M (which neglects the interaction and regenera- 
tion terms, since Am A M at low energies) and the high energy solution <p M (which 
neglects instead the decay term, since Am <§C A m at very high energies) [2, 3, 4]. 

The high energy solution is given by 

Zktm p~ x / a m _ P -X/A N 

K(E,X) = Znm ; faE-i , (11) 

where the meson attenuation length Am is defined analogously to equation (10). 

For the low energy case, the solution is obtained by replacing A M —> \ d M in 
the last equation. Recalling that A M cx E, the resulting expression couples the 
dependence on the energy and the slant depth, and hence it does not lead to just 
a simple power law energy spectrum. However, for not very small values of X the 
solution reduces to 

<fa(E, X) = -^-^e-^^E-i . (12) 

In equations (9)-(12), Znn, Znm, and Zmm are to be evaluated for the nucleon 
spectral index 7. Notice also that the high energy meson flux has the same spectral 
index as the nucleon flux, while the low energy meson solution is flatter by an extra 
power of the energy, due to the energy dependence implicit in \ d M . 

3.2 Neutrino flux 

The neutrino flux produced via the meson weak decay can be calculated following 
the same scheme described above for the atmospheric nucleon and meson fluxes. 
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Indeed, the neutrino flux is given by a transport equation analogous to equation 
(4), but actually much simpler since it contains only the source terms arising from 
meson decay, i.e. 

^(E,X)=Y,Smu(E,X), (13) 

where 

9 (F X)- r c]F > ME',X) 1 dT Mv (E, E') 

SmAE,X) - J e dE xUEf) Tm{e) . (14) 

The decay distribution can be put in terms of F Mu , the inclusive neutrino spectrum 
in the decay of meson M, by means of the relation 

TM(E) dTM dE E) = B *"> F ">>( E > ^ > ^ 15 ) 

where Bmu is the branching ratio for the decay of meson M into a state with the 
given neutrino v. In the ultrarelativistic limit, the inclusive neutrino spectrum 
scales as F Mu (E, E') = F Mv (E/E')/E' [2, 3]. Moreover, we have already derived 
the asymptotic solutions for the meson fluxes, which are of the form <p M (E,X) oc 
E~@g(X), with f3 — 7 (/3 — 7 — 1) for the high (low) energy fluxes (see equations 
(11) and (12)). Thus, the source terms given by equation (14) can be rewritten as 

Sm,(E,X) = ^^Z^(E), (16) 

where the (3— dependent meson decay moments are defined by 3 

Z P Mv {E) = B Mu f 1 dx x^Fm^x) . (17) 
Jo 

The decay moments used in this work were calculated by fitting the /3— dependence 
of the relevant moments tabulated in [4] for 2.7 < (3 < 4, which were determined 
with Lund Monte Carlo simulation programs. 

Then, the asymptotic low and high energy neutrino fluxes can now be determined 
from the equations (13) and (16), which involve the meson decay moments and the 
corresponding high and low energy meson fluxes derived in equations (11) and (12). 
The low energy solution reads 

tf(E,X) = 2T Mv -^- (l - e-*/ A ") <f> 0N E-i . (18) 

1 — ZjnN v ' 



3 Thc branching ratio is here included in the definition of the decay moments, in analogy with the 
production/regeneration moments that include the multiplicity of the final states. This coincides 
with [4] but differs from [3]. 
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As expected, the neutrino flux develops rapidly, on the scale of a nucleon interaction 
length, and then remains stable. Hence, at ground level (X S> An) the low energy 
flux will be 

tfW = ZLuT^y-faBn . (19) 
J- — Ann 

Similarly, the high energy flux at ground is 

A,H (F \ _ 77+1 Z nm ln(A M /Ajy) e M , p -( 7 +i) r9m 

where 9 is the zenith angle of the CR incidence direction considered and 

m M ch 

£m — , l^-U 

with rriM and Tm the meson's mass and mean life, respectively, and ho = 6.4 km 
a typical scale height for density variations in the atmosphere. Since the produc- 
tion/regeneration moments Z NN , Z NM , and Z MM are the same as those appearing in 
the meson fluxes, they are to be evaluated for the nucleon spectral index 7. Equation 
(20) is actually valid for zenith angles sufficiently small so that the curvature of the 
earth can be neglected, and assumes a density profile of an isothermal atmosphere 
(i.e. p{h) oc exp(—h/h )), which is appropriate in the stratosphere (h > 11 km) 
where most particle interactions occur [3, 4]. For very large angles, one can still 
use equation (20) with the simple prescription of replacing 9 — >■ 9*(9,h = 30 km), 
where for a given line of sight that corresponds to an angle 9 at the observer's 
position, the angle 9*(9,h) is the zenith angle that would be observed at a loca- 
tion where this same line of sight is at a height h with respect to the ground [3]. 
For instance, for CRs incident in the horizontal direction (9 = 90°), one has that 
9* ~ arcsin(l + h/R^)' 1 ~ 84.5°. 

In order to obtain a general expression for the atmospheric neutrino fluxes deep 
in the atmosphere, the asymptotic solutions given by equations (19) and (20) can 
be joined together by means of the interpolation function 

*-wSs' (22) 

which is actually analogous to the interpolation function given in equation (1) to 
join the fluxes associated to each galactic CR component below and above their 
knees. 

We described so far a calculational scheme that provides the nucleon, meson 
and neutrino fluxes produced by an initial power law nucleon flux <p N (E,X = 0) = 
4> 0N E"' y hitting the top of the atmosphere. We now seek the corresponding results 
for the full CR spectrum arriving at the earth, taking into account both the CR 
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composition and the break in the CR spectra at the knee. Let <f>z be the CR flux 
associated to the CR component of nuclei of charge Z and mass Am^ (with A being 
the average mass number corresponding to element Z). This nuclear component 
provides a nucleon flux 4>n,z given by <Pn,z(E n ) = A(p z (E = AE N ). If the flux 
of the CR component of charge Z is <pz{E) = <poz{E / 'E^)" lz , the corresponding 
nucleon flux is then given by 



/From these considerations, the procedure to calculate the neutrino flux produced by 
the galactic CR flux parametrized by equations (1) and (2) is now straightforward. 
Indeed, the CR component of charge Z is represented by two power law fluxes 
<p z and (f>z, with spectral indices az and az + Aa, respectively. For each power 
law flux, we can compute the corresponding nucleon flux given by equation (23) 
and use it to calculate the associated neutrino flux. Then, the final neutrino flux 
produced by the given CR component can be obtained by interpolating the two 
solutions, using again the interpolation function analogous to equation (1). Finally, 
the total neutrino flux results from summing over the contribution of all galactic 
CR components together with the additional contribution from the extragalactic 
component given by equation (3). 

Figure 2 shows the atmospheric neutrino fluxes that correspond to the Aa = 
2/3 rigidity dependent scenario, with the parametrization formulae and parameter 
values given in section 2. The figure shows the contribution coming from different 
components and for different cases, namely the vertical and horizontal conventional, 
prompt charm/GBW and prompt charm/CTEQ3 fluxes. For the prompt charm 
flux calculated with the GBW model, the figure also shows the total vertical and 
horizontal neutrino fluxes. Finally, in the latter case the conventional, charm and 
total neutrino fluxes produced by the extragalactic CR component is also indicated. 
As is clear from equation (23), the induced neutrino fluxes are sensitive, on the one 
hand, to the CR composition assumed and, on the other hand, to the model assumed 
to reproduce the steepening of the CR spectrum at the knee. Previous works that 
have taken the knee into account have assumed for simplicity that the dominant CR 
component consists of protons alone, having a spectral change Aa = 0.3 at E knee . 
In order to observe the effect of the CR composition on the neutrino fluxes, figure 3 
compares the total horizontal atmospheric fluxes (with the prompt charm component 
calculated with the GBW model) as obtained for different assumptions on the CR 
spectrum, namely the Aa = 2/3 rigidity dependent scenario (that corresponds to 
the results given already in figure 2), the same total CR spectrum but assuming 
it consists only of protons, and the CR spectrum used in [4, 7] (which is assumed 
to be constituted by protons and has a single spectral index change of Aa = 0.3 
at E knee = 5 x 10 15 eV). The latter CR spectrum has a normalization which is 




(23) 
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Figure 2: Atmospheric (i^ + P^) fluxes for the Aa = 2/3 rigidity dependent scenario 
(that corresponds to the CR spectrum plotted in figure 1). The vertical and horizon- 
tal conventional, prompt charm/GBW and prompt charm/CTEQ3 fluxes are shown, 
as well as the total vertical and horizontal neutrino fluxes for prompt charm/GBW. 
In the latter case, the conventional, charm and total neutrino fluxes produced by 
the extragalactic CR component is also indicated. 

lower by a factor of ~ 2.6 compared to the former ones, and indeed it seems to 
lay somewhat below the present observational data. Comparing the results for the 
same total CR spectrum, one can observe that the neutrino flux corresponding to 
the composition of different nuclear species is clearly below that produced by a CR 
spectrum formed only by protons, and also that this effect is more significant at 
higher energies. Indeed, equation (23) shows that the nuclear CR component of 
charge Z and mass A gives a contribution to the nucleon flux which is suppressed 
by a factor of A 2 ~ az ~ A~ ' 7 below the respective knee, while it becomes suppressed 
by a factor A 2 ~ a z-&a ^ A~ 1A we n above it. This effect was not appropriately 
treated in [6], where the suppression was considered to be given by a factor of A~ 2 , 
and moreover the impact of the change in the composition at the knee was never 
discussed before. The strong suppression of the v flux produced by the heavier 
components implies that the light (H and He) components are still responsible for 
a large fraction of the neutrinos above their knee, and hence the change in slope of 
the individual components is also reflected in the change in slope of the neutrino 
fluxes. 
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Figure 3: Comparison of the total horizontal atmospheric (v^ + v^) flux for prompt 
charm/GBW produced by CRs under different assumptions concerning spectrum 
and composition: the Aa = 2/3 rigidity dependent scenario that considers the con- 
tribution of CR components with charge 1 < Z < 28; the same total CR spectrum 
but assuming it consists only of protons; and the CR spectrum used in [4, 7] (TIG), 
which also takes only protons into account but has a lower normalization. 

4 The flux of galactic neutrinos 

The neutrino fluxes produced during the propagation of the CR particles through the 
ISM can be determined following the procedure described in the preceding section for 
the production of atmospheric neutrinos. For simplicity, we consider that the CRs 
are distributed homogeneously within the galactic disk. The highly relativistic CR 
particles interact with the ambient gas in the ISM, which is a very low-density, non- 
relativistic plasma constituted chiefly by atomic and molecular hydrogen that can be 
described by an homogeneous nucleon density ujsm — 1/cm 3 [35]. The treatment 
follows essentially the same considerations described above, but in this case one 
has to consider that the interaction of the secondaries produced in the nucleon- 
nucleon interactions is completely negligible due to the extremely low ISM densities, 
which allow for their decay well before interactions can take place, irrespective of 
the energy considered. Hence, the main contributions to the diffuse neutrino flux 
produced in the Galaxy arise from the decay of pions and muons, while the additional 
contributions coming from production and decay of heavier mesons can be safely 
disregarded. 
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As in the preceding section, let us first consider an initial nucleon flux given by 
<Pn(E,X = 0) = (PqnE" 1 . Notice that in this context the atmospheric slant depth 
should be replaced by the ISM column density traversed along the line of sight (i.e. 
X — > Rm N n IS M, where R is the distance measured from the border of the galactic 
disk along the line of sight). Due to the extreme faintness of the ISM, we should 
now take into account that 4 X <C X N . Indeed, it turns out that 

where the nucleon-nucleon total cross section is 



&NN 



35.49 + 0.307 In 2 (s/28.94 GeV 2 )] mb , (25) 



and where s ~ 2m N E stands for the center of mass energy squared [36]. 

Hence, considering equation (18) with X <C \n, the neutrino flux produced from 
pion decay in the ISM is 

ME, X) = Zl v Z NM ^4> m E-i . (26) 

an 

In order to determine the muon decay contribution, it suffices to estimate it 
directly from the results obtained already for pion decay. From the decay kinematics, 
the mean fraction of energy (relative to the parent particle) in the 7r — > n + 
decay is K w = 0.21 for the and 0.79 for the muon [1]. Analogously, in the 
/i — > e + u e + decay the effective fraction of energy for the resulting is 0.35 [1], 
i.e. a fraction = 0.28 relative to the original pion. Hence, one expects a muon 
decay contribution approximately a factor (K^/K^) 1 larger than the neutrino flux 
from pion decay (for instance, a factor of 2.1 for a nucleon spectral index 7 = 2.7, 
in sensible agreement with more detailed calculations [7]). 

As commented above, the galactic neutrino flavor oscillations may redistribute 
the produced and u e fluxes evenly among the three flavors. In that case, the 
calculation of the electron neutrino flux produced by CR interactions with the ISM 
is also required in order to determine the galactic muon neutrino flux reaching the 
earth. This flux is mainly produced in the fi — > e + v e + decay, in which the 
effective fraction of energy for the v e is 0.3 [1]. Thus, the v e flux is approximately 
a factor (0.3/0.35) 7 times the flux produced in the muon decay (i.e. a factor of 
~ 0.7 for a nucleon spectral index 7 = 2.7, again in good agreement with previous, 
more detailed calculations [7]). Hence, if the neutrino flavors are redistributed by 

4 In this case, Xn corresponds to the interaction length for nucleons propagating through the 
ISM. Notice that Aat, expressed in g/cm 2 , is actually very similar to the one corresponding to 
nucleon propagation in air. In contrast, the meson decay length Xff, which is linear in the mass 
density of the medium, will be much smaller in the ISM than in the atmosphere. 
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Figure 4: Neutrino fluxes {y^ + P M ) produced by CRs in the ISM (for different 
path lengths traversed through the ISM, namely R = 1 and 20 kpc) and in the 
atmosphere (for the horizontal direction, using the GBW model for the prompt 
charm contribution) assuming different CR compositions: (a) the Act = 2/3 rigidity 
dependent scenario that considers the contribution of CR components with 1 < Z < 
28; (b) the same total CR spectrum but assuming it consists only of protons. For 
comparison, also shown are the Waxman-Bahcall predictions for the astrophysical 
diffuse neutrino fluxes produced in their model of GRBs. 

oscillations, the (y^ + P M ) galactic flux will be a factor (1 + 0.7 x 2. 1/3. l)/3 ~ 1/2 
of that depicted in figure 4, which neglects the neutrino oscillations. 

In order to calculate the galactic neutrino flux induced by the full CR spectrum, 
the procedure follows the same steps already described above for the atmospheric 
fluxes. The total galactic neutrino fluxes are shown in figure 4(a) for R = 1 and 20 
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kpc, and correspond to the Aa = 2/3 rigidity dependent scenario (with the same 
parametrization used for the previous figures and detailed in section 2). For com- 
parison, the figure also shows the total atmospheric neutrino flux in the horizontal 
direction with the prompt charm component calculated using the GBW model, as 
well as the Waxman-Bahcall flux predictions for the astrophysical diffuse neutrino 
fluxes produced in their model of GRBs [12, 13]. Analogously, figure 4(b) shows the 
results corresponding to the same total CR spectrum but assuming it to consist only 
of protons. In agreement with previous results [7], we observe that the galactic flux 
in the direction orthogonal to the plane (corresponding for instance to R — 1 kpc) 
remains below the atmospheric flux up to the highest energies. Along directions near 
the galactic plane (e.g. for R = 20 kpc), the flux of galactic neutrinos produced in 
the ISM instead overcomes the atmospheric flux at energies larger than ~ 10 14 eV. 
The energy at which the galactic flux actually dominates depends to a large extent 
on the particular model adopted for charm production, as discussed above, and the 
results naturally also vary according to the assumed ISM column density and the 
zenith angle of arrival direction in the atmosphere. These figures show very clearly 
the significant effect of suppression in the neutrino fluxes that results from consid- 
ering a rigidity dependent scenario. Indeed, above ~ 10 15 eV the Waxman-Bahcall 
flux is found to overcome the background of galactic neutrinos with trajectories 
nearly contained in the galactic plane. However, if the same total CR spectrum 
is taken to be formed by protons alone, the resulting background flux of galactic 
neutrinos turns out to be larger than the Waxman-Bahcall flux up to ~ 10 16 eV, 
and then both fluxes stay with comparable magnitude up to the highest energies. 
Hence, rigidity dependent scenarios result in a significantly reduced background for 
the search of astrophysical neutrino sources of current interest, such as for instance 
GRBs or AGNs. 



5 Conclusions 

We computed in this work the high energy diffuse fluxes of muon neutrinos (and 
antineutrinos) produced by CRs hitting the upper atmosphere or interacting with 
the ISM in the Galaxy, showing that the results are quite sensitive to both the total 
CR spectrum and the CR composition assumed. Concerning the former, we com- 
pared previous results (which adopted an overall normalization for the CR spectrum 
actually somewhat below present observational data, and a single and abrupt slope 
change representing the knee) with those obtained from considering a smooth CR 
spectrum that soundly fits the observations and reproduces its main relevant fea- 
tures (namely, the knee, the second knee and the ankle). With respect to the latter, 
we showed that taking into account a CR composition that turns heavier above the 
knee (i.e. in agreement with the scenarios that explain the knee as due to a rigidity 
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dependent effect) the induced neutrino ffuxes become significantly suppressed, thus 
making their detection harder but also reducing the background for the search of 
other astrophysical neutrino sources. In particular, we compared our results with 
the Waxman-Bahcall predictions for the neutrino fluxes produced in their model of 
GRBs, and found indeed that the signal-to-background ratio for this model becomes 
significantly enhanced when considering a CR composition with the contribution of 
different nuclear species, as opposed to the case in which the full CR spectrum is 
assumed to be composed of protons alone. An analogous effect would lower the 
neutrino flux produced by the extragalactic component (which gives the dominant 
contribution to the total neutrino flux above E ~ 10 17 eV) if it were not constituted 
just by protons, as we assumed here. It is then clear that if the diffuse neutrino 
fluxes predicted in the GRBs models [12] or the even larger ones associated to some 
AGN models [11] do actually exist, the observation of atmospheric or galactic CR 
induced neutrino fluxes above 10 14 eV will remain probably hopeless. 

Let us also mention that it may be of interest to consider the atmospheric muon 
fluxes produced by CRs. For instance, it has recently been suggested [28] that the 
observation of the down-going atmospheric muons with neutrino telescopes would 
provide an indirect measure of the prompt atmospheric neutrino flux, and hence this 
could be used to confront the NLO QCD predictions. Indeed, due to the charmed 
particle semileptonic decay kinematics, it turns out that the prompt muon flux 
coincides with the prompt neutrino flux to within ~ 10%, irrespective of the energy 
and independently of the model used to treat the atmospheric charm production 
[28]. On the other hand, the atmospheric conventional muon flux is about a factor 
of ~ 5 larger than the conventional muon neutrino flux within the energy range 
of interest of this work, and exhibits roughly the same energy dependence [28]. 
Completely analogous effects as those described here for the neutrino fluxes should 
then be expected for the atmospheric muons produced by CRs reaching the earth. 

Finally, we stress the key importance of determining confidently the CR compo- 
sition around and above the knee, since it appears in this context as decisive in order 
to estimate reliably the diffuse high energy neutrino background produced by CRs. 
Moreover, solving the CR composition puzzle will also provide a valuable means of 
testing the different proposals concerning the origin and nature of the knee, and will 
thus shed new light on this long standing problem. 
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